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SUMMARY 


Patched-grid algorithms for the analysis of complex configurations with an im- 
plicit, upwind-biased Navier-Stokes solver have been investigated. Conservative 
and non-conservative approaches for performing zonal interpolations were im- 
plemented. The latter approach yields the most flexible technique in that it can 
handle both patched and overlaid grids. Results for a two-dimensional blunt 
body problem show that either approach yield accurate steady-state shock loca- 
tions and jump conditions. In addition, calculations of the turbulent flow through 
a hypersonic inlet on a three-zone grid show that the numerical prediction is in 
good agreement with the experimental results 1 


Through the use of a generalized coordinate transformation at the zonal interface 
between two or more blocks, the algorithm can be applied to highly stretched vis- 
cous grids and to arbitrarily-shaped zonal boundaries. Applications were made 
to the F-18 aircraft at subsonic, high-alpha conditions, in support of the NASA 
High-Alpha Research Program. The calculations were compared to ground-based 
and flight test experiments and have been used as a guide to understanding 
the ground-based tests, which are laminar and transitional, and their relation- 
ship to flight 2 3 . Calculations about a complete reconnaissance aircraft were also 
performed in order to further demonstrate the capability of the patched-grid 
algorithm 4 . 


During the period of this grant, the papers in the references section were written 
which cover in detail the three-dimensional zonal techniques that were investi- 
gated in this study. As such, they have been included as the body of this final 
report. 
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ABSTRACT 

Implicit algorithms for the compressible Navier- 
Stokes equations are applied to the computation of 
high-speed inlet flows typical of those encountered on 
propulsion systems of advanced hypersonic vehicles. 
Patched grids are used to accurately treat blunt cowl 
shock-shock interaction effects at the inlet entrance 
and efficiently compute effects due to variations in the 
downstream inlet geometry. Two algorithms, each in- 
corporating upwind-biased spatial differencing for the 
convective and pressure terms within a flux-difference- 
splitting framework, are compared to experimental 
data for internal compression ratios of two and eight 
at a nominal Mach number of 7.4. Results from the 
two codes, each with an algebraic turbulence model, 
agree closely with each other and predict the over- 
all boundary layers of velocity and temperature rea- 
sonably well, although a definitive assessment is not 
possible because of uncertainties in the experimental 
data accuracy. 

NOMENCLATURE 


A 

Jacobian matrix 

F,G 

inviscid fluxes 

F,.,G,. 

viscous fluxes 

J 

transformation Jacobian 

L 

reference length, also difference operator 

M 

Mach number 

Pt 

Prandtl number, taken as 0.72 

Q 

conservation variables 

Rei 

Reynolds number, Rei = PooU~*jL / V-yj 

T 

temperature, also diagonalizing matrices 


U, V contravariant velocities 

o speed of sound 

e total energy per unit volume 

p static pressure, p = ( 7 - l)[e - p(u 2 + v 2 )/2j 

p p pitot pressure 

q vector for state interpolations 

q Xi heat flux terms 

t time 

u,v Cartesian velocities in x and y directions 
x , y Cartesian coordinates 

At time step 

A Q change in Q over time step At 

7 ratio of specific heats, taken as 1.4 

(,rj general curvilinear coordinates 

k spatial differencing parameter 

A coefficient of bulk viscosity, A + 2p/3 = 0 

A diagonal matrix of eigenvalues 

p coefficient of molecular viscosity 

p density 

T ri rj viscous shear stress terms 

Subscripts: 

x,y denotes differentiation with respect to x,y 
t denotes total conditions 

w denotes wall value 

00 denotes conditions at freestream 

Superscripts: 

denotes quantities in generalized coordinates 
denotes dimensional value or Roe-averaging 
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INTRODUCTION 


There is an increasing effort in the development 
and validation of computational methods for high- 
speed flows because of the renewed interest in the 
hypersonic flight regime such as, for example, in the 
National Aero-Space Plane project. The propulsion 
system of advanced hypersonic vehicles will likely use 
the external vehicle contours as compression and ex- 
pansion surfaces for the inlet and nozzle, respectively, 
and thus the engine-airframe integration is a critical 
consideration. The proper prediction of mass flow into 
the inlet is dependent on the accuracy to which the 
developing boundary layers on the forebody compres- 
sion surfaces can be predicted. High heating rates 
will occur on the blunt cowl lip at the inlet entrance 
and vary considerably with Mach number and inci- 
dence because of the shock-shock interactions. The 
design of such systems relies heavily on the develop- 
ment of computational fluid dynamic codes with the 
required geometric flexibility and requisite physical 
models, since many of the high- enthalpy, high Mach 
number conditions cannot be simulated with ground- 
based facilities. 

Two recently-developed algorithms for the com- 
pressible Navier-Stokes equations are applied to two 
high-speed inlets which have features representative 
of those encountered in the propulsion systems of ad- 
vanced hypersonic inlets. The inlets were tested ex- 
perimentally in the 3.5 Foot Hypersonic Tunnel at 
NASA Ames Research Center (Gnos et al., 1973) at a 
nominal fxeestream Mach number of 7.4 and Reynolds 
number of 8.86 x 10 6 per meter. The inlets tested 
were identical on the compression surfaces ahead of 
the inlet and the internal inlet geometry was varied 
to provide internal compression ratios of two, eight, 
and twelve, leading to the designations P2, P8, and 
P12 inlets, respectively. Only the P2 and P8 inlets 
are studied here since the PI 2 inlet had fluctuations 
in surface pressures which prevented detailed experi- 
mental measurements. 

The present work is a follow-on to a previous 
study (Rudy et al., 1988) in which results from central 
and upwind difference codes were compared for the P8 
inlet with a tangent wedge geometric approximation 
to the blunt cowl lip. A patched grid formulation 
is used here to accurately treat the blunt cowl inter- 
action with the forebody compression field and effi- 
ciently treat variations in inlet geometry downstream 
of the inlet entrance. Results from the two implicit 
finite-volume algorithms, each based on an upwind- 
biased spatial differencing approach for the convective 
and pressure terms within a flux-difference-splitting 
framework, are compared with experimental results 


for both inlets. The turbulence is modeled here with 
an algebraic eddy-viscosity model (Baldwin and Lo- 
max, 1978), although alternate models could easily be 
investigated, such as in the related recent work of Ng 
et al. (1988). 

SEMI-DISCRETE FORMULATION 

The two algorithms used are both implicit 
finite- volume algorithms. The first, CFL3D (Com- 
putational Fluids Lab 3-D code), was developed by 
Thomas et al. (1985,1987) for the thin-layer Navier- 
Stokes equations. The second, USA-PG2 (Unified So- 
lution Algorithm - Perfect Gas, 2-D code), was de- 
veloped by Chakravarthy et al. (1985,1988). Both 
codes solve the time-dependent two-dimensional com- 
pressible Navier-Stokes equations expressed in strong 
conservation law form and generalized coordinates as 
below. 
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The variables £ and tj correspond to the coor- 
dinates parallel and normal to the inlet walls, respec- 
tively. The equations are nondimensionalized by p , x , 
Uoo, Aoot and L. The shear stress and heat flux terms 
are defined in tensor notation for laminar flow as: 
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where b Xl = ujr XiXj - q X{ and Sutherland’s law is 

used for the variation of molecular viscosity with 
temperature. 

The equations are represented with a semi- 
discrete finite- volume formulation. Integrating (1) 
over a cell volume bounded by surfaces of constant 
( and rf gives 

_ B. . 

( dt ], ' J ~ R, ' J 

= -(F-F v ) i+1/2J + (F-F v )i_ l/2J ( 2 ) 

- (G - Gy),-j +1 /2 + (G - G„),j_ 1 /2 

where the subscripts ij denote evaluation at cell- 
centered locations corresponding to £,,77;, respec- 
tively. A uniform spacing is taken in the computa- 
tional domain, i.e. = 1, such that the cell 

volume is associated with the discrete evaluation of 
the inverse of the transformation Jacobian. 

The term R represents the steady-state residual, 
the accuracy of which is determined by the spatial 

differencing, discussed subsequently. The equations 
are advanced implicitly in time to a steady state with 
a spatially-factored algorithm in delta form 

L i L, ] AQ=-AtR (3) 

For efficiency, first-order spatial differencing is used 
on the left-hand side operators and an approximate 
diagonal inversion is used in the £ direction, aluug 
which the grid is not highly stretched. 


SPATIAL DIFFERENCING 

The spatial differencing for the convective and 
pressure terms uses a MUSCL (Monotone Upstream- 
Centered Conservation Law) approach to evaluate the 
interface flux. The interface flux in the £ direction, for 
example, is constructed along a line of constant rjj as 


F i+i/2 = F{q ) 9 + )i+ 1/2 

= F{ <?i-i . <?<• > 9 «-n . Qi+2 » ( V£/ J ) 1+1/2) 

(4) 

where q^ + , , 2 denote state variables determined from 


upwind-biased interpolations to the left and right of 
the i+1/2 cell interface location as 

(Oi+1/2 = «+{J[(* - K )^~ + C 1 + *)£+]}« 

U + )i+l/2 = ~ K J%+ + U + K )A-]}i+l 

(5) 

where 


(A+)i = q i + 1 - q; (A_); = 

The parameter k determines the order of the scheme 
in smooth regions (k = 1/3 corresponds to the third- 
order formulation and is used herein) and slope lim- 
iting is used to maintain monotonicity across discon- 
tinuities in the flow, such as shock waves, 

A+ = max[0,min(A+sgnA--,/3A_sgnA+)]sgnA+ 
A_ = max[0, min( A_sgnA+,/3A+sgnA_)]sgnA_ 

(6) 

where 0 = (3 - *)/(l - «)• The state-variable 
interpolations are done on the vector g, which in the 
case of CFL3D is taken as the primitive variables, 
q = (p, u, v,p) T - In the case of USA-PG2, the 
characteristic variables are used, q = T~ l Q, where 
corresponds to the matrix of left eigenvectors of 
the Jacobian matrix, A = dF/dQ. 

The flux is constructed through the solution of a 
local Riemann problem at the interface using the flux- 
difference-splitting approach of Roe (1981), which can 

be written in terms of the conserved variables at the 
interface as 

^.1/2 = I [ F(Q + ) + F(Q-)-\A\ (Q+-Q-) ]. fl , 

( 7 ) 

where | A\ = T | A| T~ l is evaluated with Roe-averaged 
variables such that F+ - F~ = A(Q+ - Q~), The 
algebraic form of the equations simplifies to a con- 
siderable degree when the primitive or characteristic 
variables are used for the state- variable interpolations, 
and all choices recover identical results for first-order 
spatial differencing. The scheme has the property of 
capturing strong shocks sharply and providing min- 
imal dissipation levels in wall-bounded viscous flows 
(Van Leer, 1987). 

The shear stress and heat transfer terms are 
centrally differenced, corresponding to a second- 
order spatial discretization. The code CFL3D solves 
the thin-layer form of the equations while USA- 
PG2 solves the complete equations. Both formula- 
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tions model the turbulence using the two-layer eddy- 
viscosity model of Baldwin and Lomax (1978). 


PATCHED GRID ALGORITHM 

Zonal and/or patched grids are attractive be- 
cause they can substantially reduce the difficulties as- 
sociated with the generation of grids over complex ge- 
ometries. Additionally, they generally lead to a more 
efficient algorithm by tailoring the grids to the local 
geometry, thereby eliminating inefficient grid topolo- 
gies, and by using mathematical models in each zone 
appropriate to the local physics. The two require- 
ments placed upon a patched grid algorithm are ac- 
curacy and conservation. Rai (1985) demonstrated 
these two requirements could be met by conservatively 
interpolating the flux across a zonal boundary, so as 
to preserve the spatial flux of mass, momentum, and 
energy. 

An alternate approach is taken here. Across a 
zonal interface, conserved variables are conservatively 
interpolated from zone 2, for example, onto the grid 
locations corresponding to an analytic continuation of 
the grid points in zone 1. These interpolated points 
are used to construct the interface flux at the zone 1 
boundary through the solution of an underlying Rie- 
mann problem. The interpolated variables are conser- 
vatively constructed, much as the spatial fluxes are 
interpolated in the spatially conservative approach, 
such that they reflect the correct amount of mass, 
momentum, and energy contained in zone 2. Like- 
wise, zone 1 data is conservatively interpolated onto 
an analytic continuation of the grid points in zone 2 
in order to construct the interface flux at the zone 
2 boundary. The algorithm is constructed so that if 
a geometric match between the grid points exists at 
the zonal interface, the scheme recovers identically the 
interior scheme discussed previously. Thus, there can 
be a slight mismatch in the spatial flux across a zonal 
boundary, on the order of the local truncation error 
of the scheme. The advantage of this approach is that 
it can readily be extended to more general situations, 
such as overlapped grids where the zonal boundaries 
need not be coincident. 

The effectiveness of the approach in maintain- 
ing conservation and accuracy across the interface is 
demonstrated for a supersonic blunt body calculation 
similar to that used previously by Rai (1985). A sin- 
gle grid and a three-zone grid are shown in fig. 1. 
Pressure contours for a freestream Mach number of 2 
are shown in fig. 2 and show remarkable correspon- 
dence between the computed shock locations for both 


cases; for the 3-zone grid, a zone boundary crosses the 
shock position near the plane of symmetry and near 
the downstream boundary of the grid. 

Figure 3 shows pressure contours for two dif- 
ferent zonal boundary locations, one upstream of the 
shock location and one downstream. Other than a 
thickening of the shock pressure contours in the re- 
gions of coarser grid zones, the zonal grids show very 
little sensitivity to the zonal boundary. 

INLET DESCRIPTION 

The zonal topology used for the inlet computa- 
tions is shown in fig. 4 for both the P8 and P2 inlets. 
Zones 1, 2, and 4 are H-type grids clustered in the 
region near the wall, containing 101 x 101, 12 x 51, 
and 151 x 101 grid points, respectively. The normal 
spacing at the wall at the inlet entrance station, x = 
32 in. (81.28 cm), is .0002 in. (.005 cm), resulting in 
turbulent inner-layer similarity parameter values on 
the order of unity along the walls. Zone 3 is a C-type 
grid describing the local geometry of the blunt cowl 
Up at the inlet entrance station, spanning from x = 31 
in. to x=32 in., and containing 201 x 101 grid points. 
An enlarged view of the grid near the inlet entrance 
is shown in fig. 5, where for clarity a subset of the 
complete grid in zone 3 is shown. 

The P2 and P8 inlets represent inlet configu- 
rations typical of a hypersonic airbreathing vehicle. 
The 2-D inlets were designed to provide an internal 
compression ratio of 2 and 8, respectively. The fore- 

body wedge is a 6.5 deg. wedge intended to match 
a design Mach number of 6 at the inlet entrance un- 
der the test conditions of a freestream Mach number 
of 7.4, allowing for boundary-layer displacement ef- 
fects. The centerbody geometries for the two inlets 
are identical upstream of x = 44.25 in. (112.39 cm) 
and the cowl geometries are identical upstream of x 
= 34.25 in. (87 cm). Since the inlet entrance flow is 
predominantly an attached supersonic flow, the zonal 
topology lends itself to an efficient solution procedure. 
For example, the solution in zone 1 can be obtained 
first. Then the solutions in zones 2-3 can be obtained 
simultaneously. These zones then serve as initial value 
solutions for the solution in zone 4 for either inlet. In 
addition, the solution for the P2 inlet was obtained 
by restarting from the P8 flowfield. Although Navier- 
Stokes solutions are obtained here, a further savings 
could be effected by substituting paraboUzed approx- 
imations for the attached weak-interaction regions in 
zones 1, 2, and 4. 

The Reynolds number based on freestream con- 
ditions is 8.86 x 10 6 per meter. The total pressure is 
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4.14 x 10 6 N/m 2 , the total temperature is 811 deg. 
K, and the wall temperature is 303 deg. K. In the ex- 
periment, transition to turbulence occurred on both 
the centerbody and the cowl surfaces for the P8 in- 
let. Transition regions were specified to correspond 
approximately to the exp^'mental values. The tur- 
bulent eddy viscosity was linearly increased from zero 
to the fully- turbulent value over a length of 5.9 in. 
(15 cm) along the centerbody beginning at 10.8 in. 
(27.5 cm) and over a length of 3.94 in. (10 cm) along 
the cowl surface beginning at 39.8 in. (101 cm). Al- 
though not shown, the transitional computations and 
the experiment both indicate a boundary-layer thick- 
ness of about 0.5 in. at the inlet entrance station. 
These same transitional locations were used for the 
P2 inlet, although analysis of the velocity profiles in 
the inlet (Gnos et al., 1973) indicates that the bound- 
ary layers remain laminar on the cowl surface. All 
of the calculations were made with a 7 of 1.4, corre- 
sponding to equilibrium conditions at the tunnel test 
section. 

P8 INLET RESULTS 

Mach number contours in the region of the blunt 
cowl lip at the inlet entrance station are shown in 
fig. 6 and pressure contours at the blunt cowl lip are 
shown in fig. 7. The forebody wedge shock passes 
smoothly through the zonal interface between zones 1 
and 3. Because the grid in zone 3 is much finer than 
that in zone 1, the pressure contours coalesce slightly 
downstream of the zonal interface, owing to a finer 
resolution of the shock. The inlet forebody shock is 
much weaker than the detached shock that develops 
on the blunt cowl lip and there is consequently only a 
minor interaction on the upper surface of the cowl. 
The detached shock on the underside of the cowl 
passes smoothly across the interface between zones 
3 and 4, as do the Mach contours associated with the 
boundary layer development on the cowl. Although 
not shown, the results agree closely with Schlieren 
photographs from the experiment. 

Computed and measured pressures along the in- 
let cowl and centerbody are shown in fig. 8; den- 
sity and pressure contours of the inlet fiowfield are 
shown in fig. 9. The computed pressures agree closely 
with the experiment over the centerbody, including 
the pressure rise associated with the impingement of 
the cowl shock onto the centerbody at x=44 in. The 
computed pressures on the cowl are higher than those 
measured in the region before the impingement of the 
reflected shock from the centerbody onto the cowl sur- 
face at x=49 in. The reflected shock impingement is 


predicted slightly downstream of that measured, al- 
though the overall pressure rise through the interac- 
tion is well predicted. From previous work, it is known 
that the impingement location is very sensitive to the 
boundary layer development on the centerbody. 

Predicted pitot-pressure and total-temperature 
profiles are compared in fig. 10 with experimental 
data. In fig. 10(a), both codes predict a pitot pres- 
sure ratio of .025 in the in viscid core region above the 
boundary layer and below the cowl shock. This value 
is higher than the experimental value of .021. Ng, 
Benson, and Kunik (1986) demonstrated close agree- 
ment with the experimentally-measured level in that 
region using a 7 of 1.38. However, this choice of 7 
cannot be justified based on the temperature present 
in the experiment, as the freestream value of tem- 
perature is no more than approximately 100 deg. K. 
Thus, the ratio of specific heats would be expected to 
be 1.4 in the inviscid portions of the flow. Real-gas 
effects may be important, however, in the flow field 
induced by the presence of the pitot probe since the 
temperature increases rapidly across the shock ahead 
of the tube to a value near that of reservoir conditions. 
Pitot-tube pressures were calculated with equilibrium 
real-gas relations, but only a slight effect was found; 
the effect of the equilibrium real-gas assumption be- 
comes more important at higher stagnation tempera- 
tures. Additional supporting evidence is provided by 
the equilibrium air calculations by Ota et al. (1988) 
which indicate constant specific heats throughout the 
inlet. 

The Mach number at the cowl inlet station out- 
side of the wedge shock determined from real-gas pitot 
formulas using the measured pitot pressure data is 
716. This indicates a longitudinal variation in Mach 
number of approximately 0.01 per foot, which is in 
agreement with the wind-tunnel calibration measure- 
ments (Appendix A of Gnos et al. (1973)). Agree- 
ment of the computed and experimental pitot pres- 
sures may require incorporation of the longitudinal 
Mach-number gradient into the calculations. In addi- 
tion, the measured total temperatures do not recover 
to the freestream total temperature in the inviscid 
core regions of the inlet. 

With these discrepancies in mind, examination 
of the profiles indicates generally good agreement of 
the boundary layer thicknesses of temperature and 
pitot pressure on both the centerbody and the cowl. 
The location of the cowl shock at x=41 in. is predicted 
in better agreement with the experiment than the pre- 
vious result (Rudy et al., 1988) which used a tangent 
wedge approximation for the cowl. The two compu- 
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tations agree closely with each other and predict the 
overall trends in pitot pressure and temperature at 
x=44 in., corresponding to the cowl shock impinge- 
ment on the centerbody, and at x=47 in., just down- 
stream of the cowl shock impingement. The pitot 
profiles at x=49.5 in. reflect the fact that the com- 
puted impingement of the reflected shock from the 
centerbody onto the cowl is predicted downstream of 
the experiment, to a somewhat lesser extent with the 
computations from USA-PG2 than with CFL3D. 

P2 INLET RESULTS 

Comparisons of the computed and experimental 
results for the P2 inlet are shown in figs. 11-13. The 
centerbody pressure comparisons indicate that the lo 
cation of the pressure increase at x=44.5 in. asso- 
ciated with the cowl shock impingement is predicted 
accurately but the overall pressure rise is substantially 
higher than the experiment. The cowl pressures show 
little variation because the reflected shock from the 
centerbody passes downstream beneath the cowl sur- 
face, as can be seen from the flowfield contours in fig. 
12, which are in good agreement with Schlieren pho- 
tographs from the experiment. Also, the cowl shock 
position evident in the profiles at x=40 in. and x=43 
in. is well predicted. 

As for the P8 inlet, given the discrepancies in 
measured pitot pressures and total temperatures, the 
boundary layer thicknesses are generally predicted ac- 
curately. The decrease in the boundary layer thick- 
ness downstream of the cowl shock impingement is 
well predicted. The predicted boundary layer thick- 
ness on the cowl is slightly greater than the exper- 
imental results, associated with the observation that 
the experimental boundary layer remained laminar on 
the cowl surface. Calculations made with a laminar 
model on the cowl predicted a slightly thinner bound- 
ary layer, with a resultant improvement in the com- 
parisons. 

CONCLUDING REMARKS 

Two implicit algorithms for the compressible 
Navier-Stokes equations have been applied to the 
computation of high-speed inlet flows typical of those 
which will be encountered on propulsion systems of 
advanced hypersonic vehicles. Patched grids have 
been used to accurately treat blunt cowl shock-shock 
interaction effects at the inlet entrance and efficiently 
compute effects due to variations in the downstream 
inlet geometry. The algorithms considered are sim- 
ilar, each incorporating upwind-biased spatial differ- 
encing for the convective and pressure terms within a 


flux-difference-splitting framework. The results from 
the two codes agree closely, lending confidence that 
the equations and turbulence model are being solved 
correctly with sufficient resolution, which is a neces- 
sary step in the code validation process. The compar- 
isons with experimental data for internal compression 
ratios of two and eight at a nominal Mach number 
of 7.4 indicate that the overall features of the inlet 
flowfield are predicted reasonably well, including the 
boundary layer thicknesses of velocity and tempera- 
ture. However, because of uncertainties in the exper- 
imental data accuracy, a definitive assessment of the 
accuracy of the turbulence model is not possible. 
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zonal interface 


Figure 1. Single and 3-zone grids for blunt body 
solutions. 


Figure 2. Pressure contours for single and 3-zone grids 
for blunt body solutions; Mx> = 2.0 . 




Figure 3. Effect of variations in zonal boundary 
locations for 3-zone blunt body solutions; Moo = 2.0 



x*0. 31.9 in. 53 6 in. 

(a.) Sketch of hypersonic inlet configuration for P8 
inlet. 
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Figure 4. - Hypersonic inlet configuration and zonal 
topologies for patched grid computations. 




Figure 5. Computational grids near inlet entrance. 


zone 1 



zone 3 







zone 4 

Figure 6. 

zone 2 

Mach number contours in region 

of inlet 


entrance computed with CFL3D Navier-Stokes code. 




(b.) Cowl. 

Figure 8. Computed and measured pressures for P8 
inlet. 



(a.) Pressure contours. 



(b.) Density contours. 

Figure 9. Flowfield contours for P8 inlet computed 

with CFL3D Navier-Stokes code from x=43 in. to 
x=50.3 in. 


Figure 7. Pressure contours in the region of the blunt 
cowl lip at inlet entrance station computed with USA- 
PG2 Navier-Stokes code. 
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15 


o Experiment (Gnos et al., 1973) 

— Navier-Stokes CFL3D 

— Navier-Stokes USA-PG2 



40 42 44 46 

x, in. 

(b.) Cowl. 


48 



(a.) Pressure contours. 



(b.) Density contours. 


Figure 11. Computed and measured pressures for P2 Figure 12. Flowfield contours for P2 inlet computed 
inlet. with USA-PG2 Navier-Stokes code from x=41.6 in. 

to x=47.6 in. 



Pu (a.) x=40 inches. (101.60 cm) 

Figure 13. Computed and measured pitot and total temperature profiles for P2 inlet. 
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(b.) x=43 inches. (109.22 cm) 
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(c.) x=45 inches. (114.30 cm) 


Figure 13. Continued. 





